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We report on recent results for the pion matrix element of the twist-2 operator corresponding to the average 
momentum of non-singlet quark densities. For the first time finite volume effects of this matrix element are 
investigated and come out to be surprisingly large. We use standard Wilson and non-perturbatively improved 
clover actions in order to control better the extrapolation to the continuum hmit. Moreover, we compute, fully non- 
perturbatively, the renormalization group invariant matrix element, which allows a comparison with experimental 
results in a broad range of energy scales. Finally, we discuss the remaining uncertainties, the extrapolation to the 
chiral limit and the quenched approximation. 



1. INTRODUCTION 

Parton distribution functions of hadrons are in 
the focus of several theoretical and experimental 
investigations. In recent years lattice QCD cal- 
culations have provided estimates for the lowest 
moments of parton distribution functions, but so 
far most of the results have been obtained at one 
fixed lattice spacing only and using perturbative 
renormalization factors for the bare operators. 

The scope of our current investigation is the 
pion matrix element of the twist-2 operator cor- 
responding to the average momentum of non- 
singlet quark densities. This operator belongs to 
two irreducible representations of the lattice i?(4) 
group. We will concentrate solely on the matrix 
element of the operator 
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since it can be computed at zero external pion 
momentum and thus provides a reliable signal. 
In particular we use Schrodinger functional (SF) 
boundary conditions ^ for our computation, 
which allows the extraction of the matrix ele- 
ment from one large plateau around the middle 
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of the time extent T of the lattice. The corre- 
lation function of the matrix element /m is ob- 
tained by inserting the O44 operator between two 
SF pion states Oq and Ot at the time boundaries 
t ^ and T and suitable normalization with the 
boundary-to-boundary correlator /i |2l3j 
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The bare pion matrix element is then defined by 
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l^^yare ^ ^ —fM{t)/h , 

considering a plateau region around T/2 where 
the relative excited state contribution e is below 
0.4%. The pion mass is obtained from the time 
dependence of the improved axial- vector correla- 
tor following [2] and requiring e smaller than 
0.1%. 

2. FINITE VOLUME EFFECTS 

Finite volume effects for pion and nucleon ma- 
trix elements have not been studied before, but 
might influence the precise determination of mo- 
ments of parton distribution functions from lat- 
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Figure 1. Finite volume effects for mj^a 



Figure 2. Finite volume effects for (x) 



tice QCD. We have investigated the effects of lim- 
ited lattice extent in detail at (3=6.1 (a=0.079 
fm) with the non-perturbatively improved clover 
fermion action and lattice sizes ranging from 
123 X 36 to 323 X 56. 

While the pion mass reaches a constant value 
already for L ~ 1.2 fm and T ~ 2.8 fm (see fig- 
ure ^ , the bare pion matrix element in figure |21 
shows very strong finite volume effects. For the 
largest quark mass L ~ 1.6 fm and T ~ 3.2 fm 
seems to be sufficient, whereas an even larger spa- 
tial extent of 1.9 fm is needed for the smallest 
quark mass. The lines in both figures represent 
exponential fits with the purely phenomenologi- 
cal ansatz X{L)=co + ci/L^^^ exp{—C2L), where 
X = niTra or (x). Please note that a power-law 
fit ansatz would describe the data almost equally 
well. 

These observations lead to the expectation that 
the finite volume effects for nucleon matrix ele- 
ments could be even larger, since the influcence 
of the limited lattice size on the nucleon mass has 
already been observed to be much stronger than 
for the pion mass 21 . This might be one ingredi- 
ent for the deviation of current lattice QCD esti- 
mates of {x) for the nucleon in comparison with 
experimental data. 

3. CONTINUUM EXTRAPOLATION 

In order to improve the continuum extrapola- 
tion we use both the unimproved Wilson and non- 
perturbatively improved clover fermion action. 



which can have different 0{a) effects due to the 
usage of the unimproved O44 operator, but should 
agree in the continuum limit. We have comple- 
mented a previous data set [3] with simulations at 
/3=6.1 (243 X 42) and /3=6.45 (32^ x 72). The bare 
matrix elements were renormalized with the ap- 
propriate Z factors computed non-perturbatively 
at the same low energy scale /xq in the SF scheme 
[S]. In figure 01 we show the combined continuum 
extrapolation (4 smaller lattice spacings) of the 
renormalized matrix element 
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where the chiral extrapolation had been per- 
formed linearly in the quark mass. This yields 
a value of {xYgp{^Jio) = 0.870(33) for the renor- 
malized pion matrix element in the SF scheme. 

4. RGI MATRIX ELEMENT 

The phenomenological analysis of the experi- 
mental data is usually done in the MS scheme. 
In order to translate our fully non-perturbatively 
computed matrix element {xYgp{^o) renormal- 
ized in the SF scheme to the MS scheme, we 
need first to calculate the universal renormaliza- 
tion group invariant (RGI) matrix element. This 
is done following the complete non-perturbative 
evolution |S] in the SF scheme until it is possi- 
ble to make (at small volumes) a matching with 
perturbation theory. Using the so-called UV- 
invariant step scaling function it is then possible 
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Figure 3. Combined continuum extrapolation for 
different discretizations of the fermion action 

to eliminate any reference to the scale /io and ob- 
tain finally the RGI matrix element 

i^f"' = (^)5F (mo) ^Fw(a^o) = 0.f92(ll) . 

The RGI matrix element allows a simple conver- 
sion to any desired scheme (e.g. MS at 2GeV) 
requiring only the knowledge of the (3- and 7- 
function up to a certain order in perturbation 
theory in this scheme The resulting num- 
ber (a;)jg^(^=2GeV)=0. 265(15) has decreased 
slightly compared to the previous lattice com- 
putation 6 (x);jg5r(Ai=2.4GeV)=0.273(12), which 
has been performed only at one lattice spacing 
and with perturbative renormalization factors. 

5. LIMITATIONS 

Our new result is still higher, but already more 
precise than the latest NLO analyses of Drell-Yan 
and prompt photon ttN data j8l9| . which yield a 
combined experimental value of (x)j;jg{fj,—2GeV) 
= 0.21(2). The quenched approximation is cer- 
tainly one limitation that has to be explored in 
the future. 

Another aspect that has been investigated re- 
cently is the chiral extrapolation. It has been 
shown U] that a non-linear fit of the previous data 
[HI leads to {x)jYg — 0.24(1)(2), which is much 
more compatible with the experimental number. 
The first error is the statistical and the second 
describes the variation with the fit parameters. 



Figure 4. Comparison of chiral expolations: lin- 
ear (solid line) and non-linear X| with A=0.4, 0.7, 
1.0 GeV (dotted lines) 

In figure^ we show our results performing first 
the combined continuum extrapolation of Wilson 
and clover data at fixed mj. A linear chiral fit 
yields a consistent result {x)jjg = 0.263(14) with 
interchanged limits, while a non-linear fit simi- 
lar to 13 yields {x)jys = 0.22(1)(2) and is thus 
compatible within error with the experimental 
number. Despite these results the correct chiral 
extrapolation needs further investigations, espe- 
cially in the region of smaller and realistic quark 
masses. 
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